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ORIGIN AND DYNAMICAL SUPPORT OF IONIZED GAS IN GALAXY BULGES 
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ABSTRACT 

We combine ionized gas ([N II] A6583) and stellar central velocity dispersions for a sample of 345 galaxies, 
with and without active galactic nuclei (AGNs), to study the dynamical state of the nuclear gas and its physical 
origin. The gas dispersions strongly correlate with the stellar dispersions over the velocity range of a « 30-350 
km s -1 , such that a^/cr* « 0.6- 1.4, with an average value of 0.80. These results are independent of Hubble type 
(for galaxies from E to Sbc), presence or absence of a bar, or local galaxy environment. For galaxies of type Sc 
and later and that have cr„ < 40 km s -1 , the gas seems to have a minimum threshold of a g » 30 km s -1 , such that 
erg/er* always exceeds 1 . Within the sample of AGNs, crg/cr* increases with nuclear luminosity or Eddington ratio, 
a possible manifestation of AGN feedback associated with accretion disk winds or outflows. This extra source of 
nongravitational line broadening should be removed when trying to use a g to estimate rr* . We show that the mass 
budget of the narrow-line region can be accounted for by mass loss from evolved stars. The kinematics of the 
gas, dominated by random motions, largely reflect the velocity field of the hot gas in the bulge. Lastly, we offer a 
simple explanation for the correlation between line width and line luminosity observed in the narrow-line region 
of AGNs. 

Subject headings: galaxies: active — galaxies: bulges — galaxies: ISM — galaxies: kinematics and dynamics 
— galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 

Warm (~ 10 4 K), ionized gas is pervasive in the central re- 
gions of nearby galaxies of all morphological types. Easily de- 
tected through optical emission lines, either by imaging or spec- 
troscopy, this component of the interstellar medium has served 
as a diagnostic of the physical conditions in galactic nuclei, 
including their source of excitation and chemical abundances. 
The optical line emission also provides a powerful tracer of the 
nuclear kinematics of the galaxy, in efforts to probe the mass 
distribution of the bulge, to detect central dark massive ob- 
jects, and to study substructure such as kinematically decoupled 
cores. 

What is the dynamical state of the nuclear gas? How does it 
relate to the dynamics of the stars? Is the gas primarily rotation- 
ally supported, or do random motions also contribute? While 
these issues have been addressed in the past, they have been 
based either on investigations of individual objects or small, re- 
stricted samples. A number of authors have noted, for example, 
that in the bulges of disk (lenticular and spiral) galaxies, the 
gas often rotates slower than the stars (Fillmore et al. 1986; 
Kent 1988; Kormendy & Westpfahl 1989; Bertola et al. 1995; 
Cinzano et al. 1999; Pignatelli et al. 2001). The same phe- 
nomenon has been found in some elliptical galaxies (e.g., Cald- 
well 1984; Caldwell et al. 1986). The physical origin for this 
effect, however, is not entirely clear. An interesting possibil- 
ity is that the gas experiences nonnegligible pressure support 
(e.g., Kent 1988; Cinzano & van der Marel 1994; Bertola et 
al. 1995; Fisher 1997), an idea consistent with the detection 
of large central line widths in many bulge-dominated galaxies 
(e.g., Demoulin-Ulrich et al. 1984; Phillips et al. 1986; Bertola 
et al. 1995). Vega Beltran et al. (2001) compiled gas and stel- 
lar velocity measurements for a sample of ^40 disk galaxies to 
examine trends with Hubble type. They conclude that, with a 
few exceptions, the gas is dynamically cold, confined to a disk 
supported mainly by rotation. 



What is the role of nongravitational forces on the gas dy- 
namics, especially in objects harboring active galactic nuclei 
(AGNs)? This issue has been investigated most thoroughly in 
the context of the narrow-line region (NLR) in radio (Heckman 
et al. 1985; Smith et al. 1990; Baum & McCarthy 2000; Noel- 
Storr et al. 2007) and luminous Seyfert (e.g., Whittle 1992b, 
1992c; Nelson & Whittle 1996) galaxies. The overall consen- 
sus from these studies is that the velocity field of the NLR gas 
largely traces the gravitational potential of the bulge of the host 
galaxy. Super-virial motions can be found, but only in the mi- 
nority of objects exhibiting powerful, jetlike radio structures, 
which presumably impart additional mechanical acceleration to 
the gas. Comparing the widths of the optical nebular lines and 
the stellar velocity dispersions of a large sample of relatively 
luminous Seyfert 2 galaxies, Greene & Ho (2005a) reinforced 
the notion that the stellar potential principally governs the kine- 
matics of the NLR, with important exceptions arising in cases 
where nuclear activity is especially strong. To date, much less 
attention has been paid to more garden-variety AGNs such as 
low-luminosity Seyferts and LINERs, or to nuclei powered by 
massive stars. 

The Palomar spectroscopic survey of nearby galaxies (Ho et 
al. 1995, 1997a, and references therein) offers an excellent 
opportunity to examine some of the issues mentioned above. 
In addition to the existing extensive database of emission-line 
widths, a catalog of central stellar velocity dispersions is now 
available (Ho et al. 2009). This allows a detailed comparison of 
the central velocity dispersions of the gas and stars for a large, 
well-defined sample of galaxies spanning a wide range of Hub- 
ble types and nuclear spectral classes. This is the goal of this 
paper. 
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FIG. 1 . — Distribution of gas (a g ) versus stellar (cr,) velocity dispersion as a function of Hubble type. Unbarred and barred galaxies are plotted as filled and 
encircled points, respectively. The dashed line denotes cr g = cr, . A typical error bar is given in the lower-right corner of the first panel. 



2. SAMPLE AND DATA 



Our analysis draws from the nearly complete, magnitude- 
limited sample of 486 bright galaxies spectroscopically sur- 
veyed using the Palomar 5-meter telescope. We focus on the 
galaxies that have well-measured widths (FWHM) for [N II] 
A6583, which, as explained in Ho et al. (1997a), is the line 
of choice in the Palomar survey for profile measurements. The 
majority of the Palomar galaxies (428/486 or 88%) now have 
central stellar velocity dispersion determined directly from the 
original survey data (Ho et al. 2009). Although many of the 
galaxies have preexisting velocity dispersion measurements in 
the literature, the Palomar data constitute a self-consistent, ho- 
mogeneous set of measurements. A significant fraction of the 
galaxies (30%) have dispersions measured for the first time. Of 
the 407 galaxies that have FWHM([N II]) measurements, 372 
(91%) have stellar velocity dispersions. If we confine our at- 
tention to the sources that have FWHM([N II]) measurements 
with a quality rating of "b" or better (see Ho et al. 1997a), ex- 
cluding a handful objects with only upper limits on their line 
widths, 345 have stellar velocity dispersions. 

The above subset of 345 objects, the focus of this paper, rep- 



resents fairly the parent population of emission-line nuclei in 
the Palomar survey. The 25 E, 70 SO-SO/a, 149 Sa-Sbc, 62 Sc, 
and 30 Scd-Sm galaxies constitute, respectively, 81%, 93%, 
96%, 84%, and 48% of the parent population of emission-line 
objects. In terms of the nuclear spectral classes, the coverage 
is excellent for the AGNs: the present subsample contains 94% 
of the Seyferts, 93% of the LINERs, and 92% of the transition 
nuclefl in the parent survey. The high completeness fractions 
simply reflect the high AGN detection rate in bulge-dominated 
galaxies (Ho et al. 1997b), for which velocity dispersions are 
easier to obtain. By comparison, H II nuclei are 72% complete 
in the current sample — still a very high fraction. 

3. THE CORRELATION BETWEEN a g AND cr* 

3.1. Dependence on Hubble Type 

Figure 1 shows the overall correlation between gas (a g ) and 
stellar (er») velocity dispersions, first for the entire sample, and 
then individually for several Hubble type bins. Most of the 
narrow emission lines, especially near the core, do not devi- 
ate strongly from a simple Gaussian (see Ho et al. 1997d), and 
we define a g = FWHM([N II])/2.35. Barred galaxies, taken to 
be those designated as "SB" or "SAB" in the Third Reference 
Catalogue of Bright Galaxies (de Vaucouleurs et al. 1991), are 



See Ho et al. (1997a) for definition of the spectral classifications. 
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FIG. 2. — Distribution of gas (<r g ) versus stellar (<r„) velocity dispersion as afunction of nuclear spectral class. The symbols are coded according to Ha luminosity. 
Objects containing linear radio sources with 6 cm powers greater than 10 21 ' 7 W Hz~' are marked with boxes. The dashed line denotes a g = cr t . A typical error bar 
is given in the lower-right corner of the first panel. 



marked for emphasis. It is clear that <j g correlates strongly with 
<7*, albeit with considerable scatter. For the sample as a whole, 
the Kendall's r correlation coefficient (Isobe et al. 1986) is 
r = 1.19, which rejects the null hypothesis of no correlation 
with a probability of P nu n < 10~ 4 . Both a g and cr* cover ap- 
proximately the same rang^l of velocities, from ^30 to 350 
km s" 1 . Although the correlation is most evident for the en- 
tire sample, the same trend replicates for the various Hubble 
type bins. It is least well-defined for the elliptical and extreme 
late-type (Scd-Sm and 10) galaxies because of the small sam- 

2 The lower limit of the velocities, however, is strongly influenced by the spectral resolution of the Palomar survey, which is a K, 40 km s -1 in the red and ~120 
km s _1 in the blue (Ho et al. 1997a). Published stellar velocity dispersions often are also not very reliable for low-dispersion galaxies. Thus, we consider the region 
defined by a g < 40 km s -1 and <r* < 40 km s _1 to be quite uncertain. 



pies and limited dynamic range in velocity dispersions. Thus, 
to first order, the gas velocities seem to track the random ve- 
locities of the stars. Quantitatively, however, the gas disper- 
sions are systematically slightly lower than the stellar disper- 
sions, on average by a constant offset of ^20%. Defining 
the residuals from the a g = c* line by Act = log a g - log cr* , 
(Act) = -0.099 ± 0.010. Barred and unbarred galaxies behave 
the same, with (Act) = -0.094 ± 0.014 and -0.104 ± 0.014, re- 
spectively. At the lowest velocities, cr* <J 40 km s , nearly all 
the objects have Act > 0. This subgroup of high-Aer objects 
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FIG. 3. — Residuals of the cr g -a t correlation, Act = log a g — log <r» , as 
a function of Ha luminosity L\\ a , for (a) all AGNs (Seyferts, LINERs, and 
transition objects), including the sample of luminous Seyferts from Nelson & 
Whittle (1995), and (b) H II nuclei. Objects with linear radio sources and 
fijem > 10 21 7 W Hz~' are marked with boxes. The dashed line denotes a g = 
a - *. The solid line in panel (a) shows the best-fit regression line. A typical 
error bar is given in the lower-right corner of the panel (b). 

cannot be attributed entirely to resolution effects because only 
a small number of [N II] profiles are unresolved and including 
them does not alter the distribution of data points. Rather, it 
appears that galaxy centers reach a "floor" at a g m 30 km s" 1 . 

In detail, the different Hubble types also have similar dis- 
tributions: (Act) =-0.15 ±0.043, -0.105 ±0.018, -0.113 ± 
0.014, and -0.083 ±0.026 for E, SO-SO/a, Sa-Sbc, and Sc, 
respectively. Ellipticals appear to have somewhat larger off- 
sets than the disk galaxies, but the means of the different sam- 
ples are not statistically different according to the Student's t 
test. By contrast, very late-type objects, in our sample dom- 
inated by Scd-Sm galaxies, clearly stand apart with (Act) = 
-0.036 ±0.041, differing mildly from the other groups at a sig- 
nificance of 94%. 

3.2. Dependence on Local Environment 

Ho et al. (1997a) give two parameters to gauge the tidal in- 
fluence of nearby neighbors: (1) p ga i, the density of galaxies 
brighter than Mb = -16 mag in the object's local vicinity, and 
(2) p , the projected angular separation to the nearest neighbor 
within a magnitude difference of ± 1 .5 mag and a velocity dif- 
ference of ±500 km s" 1 . We find that Act shows no correlation 
with either p ga i or 6 p . 

3.3. Dependence on AGN Properties 

We next evaluate the dependence of the <7 g -er* correlation 
on nuclear spectral classification (Fig. 2). The four classes 

3 Whittle (1992b, 1992c) assumes a Hubble constant of H = 50 km s^Mpc -1 
reference wavelength of 6 cm. The corresponding threshold, for a spectrum /„ oc 



of emission-line objects appear broadly similar, but in detail 
there are some subtle, important differences. Taken at face 
value, the Seyfert sample shows the least degree of offset 
from the a g = c* line. Its distribution of Act, with an av- 
erage value of -0.054 ± 0.025, statistically differs from that 
of transition nuclei ((Act) = -0.13 ± 0.021) at a significance 
level of 98% according to the Student's t test. LINERs, with 
(Act) = -0.095 ± 0.019, lie sandwiched in between, although 
formally the differences are only statistically marginal. Com- 
bining all three AGN subclasses, (Act) = -0.096 ± 0.012. 

As a group, H II nuclei have an overall (Act) = -0.103 ± 
0.017. However, as mentioned in § 3.1, late-type spirals, which 
largely dominate the H II class, seem to comprise two distinct 
populations divided roughly at ct* w 40 km s" 1 . Below this 
value, our sample of H II nuclei has (Act) = +0.28 ±0.055; 
above it, (Act) =-0.15 ±0.013. 

Whittle (1992b, 1992c) and Nelson & Whittle (1996) sug- 
gested that Seyfert galaxies with radio sources that are both 
strong and linearly extended tend to have emission lines that are 
systematically broader than the stellar velocities. They interpret 
this to mean that the outflowing radio jets provide a secondary, 
nongravitational source of acceleration to the NLR gas. Whittle 
(1992b, 1992c) chose P 2 0cm = 10 22 - 5 W Hz" 1 as the threshold be- 
tween strong and weak radio sources. We have adopted this cri- 
teriorQ, and with the radio data for the Palomar low-luminosity 
AGNs presented in Ho & Ulvestad (2001) and elsewhere (Ho 
2008, and references therein), highlighted in Figure 2 the ob- 
jects with luminous, jetlike radio sources. Two of the objects in 
this subset — NGC 1068 and NGC 3079 — are clearly systemati- 
cally offset above the rest. But another six that satisfy the same 
radio criteria do not look extraordinary whatsoever, and not all 
high-Ac objects have strong radio jets. The same trends hold 
for luminous radio sources among the LINERs and transition 
objects. As a whole, strong radio jets are neither necessary nor 
sufficient to produce super-virial velocities in the NLR. 

Closer inspection of Figure 2 reveals a more interesting effect 
within the AGN samples. When the objects are flagged accord- 
ing to Ha luminosity, those with high luminosities — here, for 
illustrative purposes, chosen somewhat arbitrarily to be those 
above Lu a = 10 40 erg s" 1 — tend to lie systematically displaced 
toward higher Act compared to objects with lower luminosities. 
(We use the narrow component of the Ha luminosity, corrected 
for Galactic and internal extinction.) This is most clearly seen 
for the Seyfert sample. It is less obvious for the LINERs and 
transition objects because of the paucity of luminous sources in 
these classes (Ho et al. 2003), but the same trend is noticeable 
nonetheless. Within the Seyfert sample, (Act) =+0.004 ±0.043 
for objects with L Ha > 10 40 erg s" 1 and (Act) = -0. 10 ± 0.026 
for objects with Lu a < 10 4() erg s" 1 . The difference between 
the two means is significant at the level of 96%. 

To explore the luminosity effect further, we augmented the 
Palomar sample with 52 Seyfert nuclei from the work of Nelson 
& Whittle (1995) that were not already in the Palomar survey. 
These additional objects provide better coverage at the high- 
luminosity end. Ha luminosities were converted from H(3 lumi- 
nosities given in Nelson & Whittle (1995) and Whittle (1992a), 
assuming Ha/H/3 = 3.1 (Gaskell & Ferland 1984) and adjusted 
to a Hubble constant of Ho = 75 km s _1 Mpc _1 . Figure 3a plots 
Act against Lh« for the combined sample of AGNs, again flag- 
ging the luminous, extended radio sources. Note the nonzero 



and a reference wavelength of 20 cm, whereas we use Ho - 
v-° n , is P 6cm = 10 21 ' 7 W Hz" 1 . 



75 km s Mpc and a 
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FIG. 4. — Residuals of the Cg-cr* correlation for the AGNs (Seyferts, LIN- 
ERs, and transition objects in the Palomar survey, plus the luminous Seyferts 
from Nelson & Whittle 1995), Ac, as a function of Pg cm , the 6 cm power The 
symbols for the different sources are given in the figure legend. Objects with 
linear radio sources and P6 cm > lO 21 " W Hz" 1 are marked with boxes. The 
dashed line denotes a g = a„ . The solid line gives the best-fit regression line. 
A typical error bar is shown in the lower-right comer. 

slope. A correlation can be seen for the Seyferts alone, as well 
as for the combined sample of LINERs and transition nuclei. 
Treating log L H a as the independent variable, an ordinary least- 
squares regression fit to all the points gives 

Act = (0.047 ± 0.015) logL Hct -(1.95 ± 0.59), (1) 

with an rms scatter of 0.17 dex. The scatter reduces to 0.16 
dex with the luminous radio sources excluded. The correla- 
tion has high statistical significance: the generalized Kendall's 
t test yields r = 0.47 and P nuU < 10" 4 . Omitting the flagged 
radio sources has little effect on the results. By contrast, the 
sample of H II nuclei (Fig. 3b) does not show this effect, even 
though the H II nuclei and AGNs span essentially the same 
range in Ha luminosities. The different behavior of H II nu- 
clei compared to AGNs cannot be attributed to differences in 
morphological types. Roughly half of the H II nuclei are hosted 
in bulge-dominated (SO-Sbc; T = -2 to 4) galaxies, and these 
follow the same pattern as the rest of the sample (Fig. 3b). 

In Figure 4, we substitute the Ha luminosity with the 6 cm 
radio power. Radio data for the Palomar Seyferts come from 
Ho & Ulvestad (2001), and those for the supplementary sample 
of luminous objects come from Nelson & Whittle (1995) and 
Whittle (1992a). Data for the LINERs and transition objects are 
taken from a variety of sources, as summarized in Ho (2008). 
Again, we see a positive, significant trend, although formally it 
is weaker than that for Lhq- In this case the Kendall's r cor- 
relation coefficient is r = 0.29, with P nu u = 6 x 10~ 4 . A linear 
regression fit, calculated using the method of Schmitt (1985), 
which properly treats the few censored radio data points, gives 

Act = (0.013 ± 0.0085) logP 6cm - (0.34 ± 0.17). (2) 

The tendency for Act to increase with optical and radio power 
suggests that the AGN luminosity, over a wide range of values, 
injects an additional, systematic contribution to the line broad- 
ening. 

Lastly, Figure 5 illustrates that an even stronger correlation 
exists between Act and the Eddington ratio, which is propor- 
tional to the mass accretion rate. As in Greene & Ho (2007), 
we estimate the bolometric luminosity from the Ha luminosity, 
converting L H a to L5100 using the line-continuum correlation of 
Greene & Ho (2005b) and then adopting L bo i = 9L5ioo- 
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FIG. 5. — Residuals of the <Tg-<T* correlation for the AGNs (Seyferts, LIN- 
ERs, and transition objects in the Palomar survey, plus the luminous Seyferts 
from Nelson & Whittle 1995), Act, as a function of Eddington ratio, Lbol'^Edd- 
The symbols for the different sources are given in the figure legend. Objects 
with linear radio sources andP 6cm > 10 21 - 7 WHz -1 are marked with boxes. 
The dashed line denotes <j g = <t„ . The solid line gives the best-fit regression 
line; the dotted line shows the relation found by Greene & Ho (2005a; Equa- 
tion 4). A typical error bar is shown in the upper-right corner. 

We use the Mbh - a* relation of Tremaine et al. (2002) to ob- 
tain the black hole mass and hence the Eddington luminosity, 
L Edd = 1-26 x 10 38 (M BH /M Q ) erg s" 1 . The Kendall's r corre- 
lation coefficient now rises to r = 0.74, with P nu n < 10~ 4 . An 
ordinary least-squares regression fit with log L bo \/L Edd as the 
independent variable gives 

Act = (0.087 ± 0.01 l)logL bol /L Edd + (0.23 ± 0.042). (3) 

The correlation is surprisingly tight, with an rms scatter of only 
0.15 dex. If we omit the sources with luminous radio jets, the 
scatter reduces even further to 0.14 dex. Act seems to corre- 
late somewhat better with Lboi/^Edd than with Ln a , although 
formally, according to the / and F tests, neither the mean nor 
the variance of the residuals about the best-fit lines differs at a 
statistically significant level. 

3.4. Geometry and Kinematics of the Ionized Gas 

The gas line widths presented in this study were extracted 
from a single 2"x4" aperture centered on the nucleus, which 
corresponds to a region 200 pc x 400 pc for a typical distance of 
20 Mpc (Ho et al. 1997a). Without spatially resolved data, it is 
impossible to know a priori in any individual object whether the 
line width represents true velocity dispersion or instead arises 
from spatial smearing of a steep inner rotational gradient. Nev- 
ertheless, we can appeal to general arguments to test the sce- 
nario of rotational broadening. If the nuclear gas is confined to 
a rotating disk aligned with the large-scale galactic plane, one 
might naively expect the line width to correlate with the galaxy 
inclination angle. This effect is not seen in the Palomar survey 
(Ho 1996; Ho et al. 2003). However, as discussed by Whittle 
(1992b), this simple test is likely to be inconclusive because the 
line widths are primarily governed by virial velocities and AGN 
properties rather than by projection effects. 

Instead, Whittle (1992b), as did Heckman et al. (1989), ex- 
amined the correlation between inclination angle and the ratio 
AVrot/FWHM, where AV rot is the observed rotation amplitude 
of the large-scale disk and FWHM is the full-width at half max- 
imum of the central line profile. If the central nebular emission 
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FIG. 6. — Distribution of the ratio AV rot /FWHM([N II]) versus sin i as a 
function of Hubble type. The dashed line denotes an idealized model wherein 
FWHM([N II]) arises entirely from rotation in the galactic plane. Actual fits 
to the data, after excluding the luminous, linear radio sources (marked with 
squares), are given as solid lines; the regressions are calculated using an ordi- 
nary least-squares fit with log sin i as the independent variable. 

arises from a coplanar nuclear disk, FWHM, like AV rot , is a pro- 
jected quantity, and AV rot /FWHM will be independent of incli- 
nation. In the extreme alternative scenario where FWHM repre- 
sents purely random velocity dispersion, AVWFWHM should 
correlate positively with inclination. The results of this exercise 
are shown in Figure 6 for all the non-elliptical galaxies, using 
data for AV rot (derived from integrated H I line widths) and in- 
clination angle i from Ho et al. (1997a). In each panel, the 
coplanar rotation model is marked with the horizontal dashed 
line. The actual regression fit to the data (after excluding the 
flagged radio sources) is shown as a solid line. With the excep- 
tion of the small number of very late-type (Scd-Sm) spirals, the 
vast majority of the sample strongly disagree with the coplanar 
rotation model. 

As a final check, we also plotted the residuals Act from Fig- 
ure 1 against sin i. We find no correlation, as to be expected if 
<j g and <7* are both inclination-independent. 



4. DISCUSSION 

4. 1 . Estimating cr* from a g 

We combine new central stellar velocity dispersions (Ho et 
al. 2009) with published emission-line width measurements 
(Ho et al. 1997a) to examine the relationship between the kine- 
matics of the ionized gas and stars in the central regions of 
galaxies. Unlike most studies that primarily focus on the [O III] 
A5007 line, ours uses [N II] A6583, which offers at least two 
major advantages. [N II] has a lower ionization potential and a 
lower critical density than [O III], making it a better tracer over 
a larger radial extent of the bulge. [N II] also provides a more re- 
liable probe of the gravitational potential because it is less sus- 
ceptible to contamination by outflows or other radial motions. 



Like most narrow emission lines in AGNs, [N II] does exhibit 
line asymmetries indicative of radial flows (Ho 1996; Ho et al. 
2003), but they seem to be less prevalent than in high-ionization 
transitions such as [O III] (Heckman et al. 1981; Vrtilek & Car- 
leton 1985; Whittle 1985). Studies of AGN line profiles that in- 
clude transitions from a wide range of ionization states (Busko 
& Steiner 1992; Greene & Ho 2005a) find that [S II] AA6716, 
673 1 generally tend to be less asymmetric than [O III] (but see 
Rice et al. 2006). Since the profile of [S II] empirically matches 
well the profile of [N II] (e.g., Ho et al. 1997d), it follows that 
[N II] is also less asymmetric than [O III]. Moreover, the aver- 
age line centroid of [N II] is consistent with the systemic ve- 
locity of the host galaxy, whereas high-ionization lines such as 
[O III] are often blueshifted (Boroson 2005). Previous work that 
makes use of [N II] to probe the kinematics of the NLR include 
those of Phillips et al. (1986), Verdoes Kleijn et al. (2006), 
Zhou et al. (2006), Chen et al. (2008), and Walsh et al. (2008). 

Our principal finding is that, within the central few hundred 
parsecs of galaxy bulges, the velocity dispersion of the ionized 
gas is comparable to, but typically somewhat lower than, the 
velocity dispersion of the stars. The rough correspondence be- 
tween o g and 0% indicates that the NLR is in dynamical equi- 
librium with the stellar potential. Defining Act = log a g - log ct* , 
Act ranges from -0.15 to +0.08, with an average value of 
-0.099, or (c^/er*) = 0.80. In detail, Act shows, at best, a 
mild variation with galaxy morphology but no discernible de- 
pendence on the presence of a large-scale bar or local galaxy 
environment. By contrast, the more luminous AGNs studied 
by Whittle (1992b) seem to exhibit broader lines in barred and 
interacting galaxies. 

As an aside, we note that the late-type spirals in our sam- 
ple show an intriguing property: below 0% « 40 km s -1 , Act is 
always positive. It seems that in these late-type galaxies with 
small bulges (possibly all "pseudo-bulges"; Kormendy & Ken- 
nicutt 2004), the velocity dispersion of the warm (~ 10 4 K) 
ionized gas has a floor of <~ 30 km s" 1 . This might represent a 
minimum threshold of interstellar turbulence in the central re- 
gions of present-day disk-dominated galaxies, maintained, per- 
haps, through stellar feedback associated with their modest star 
formation rates (Ho et al. 1997c). The dispersion threshold 
we find is qualitatively similar to the magnitude of non-ordered 
motions detected in the central regions of low-surface bright- 
ness galaxies (e.g., Kuzio de Naray et al. 2008; Pizzella et al. 
2008a, 2008b). 

The availability of a uniform set of nuclear spectroscopic 
classifications for the Palomar galaxies has prompted us to re- 
examine the kinematics of the NLR, their relationship with the 
stellar potential of the bulge of the host galaxy, and the possi- 
ble role of nongravitational perturbations by AGNs. Previous 
treatment of this subject has concentrated almost exclusively 
on radio galaxies (Smith et al. 1990; Verdoes Kleijn et al. 
2006; Noel-Storr et al. 2007) and relatively luminous Seyfert 
nuclei (Veron 1981; Terlevich et al. 1 990; Veilleux 1 99 1 ; Whit- 
tle 1992b, 1992c; Nelson & Whittle 1996; Jimenez-Benito et 
al. 2000; Botte et al. 2005; Greene & Ho 2005a; Bian et 
al. 2006; Zhou et al. 2006; Chen et al. 2008; Dasyra et al. 
2008). The reviews by Wilson & Heckman (1985) and Whittle 
(1993) included some preliminary results on LINERs, although 
they were based on rather fragmentary data. They concluded 
that the NLR kinematics of LINERs, as in Seyferts, are largely 
virial. In their survey of ionized gas in E and SO galaxies, most 
of which contain LINER nuclei, Phillips et al. (1986) compared 
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the FWHM of [N II] A6583 with the absolute magnitude of the 
host galaxies and inferred that in general a g <; 0% . 

Our analysis strongly reinforces the notion that gravity dom- 
inates the kinematics of the NLR, not only in Seyferts, but in- 
deed in all emission-line nuclei found in nearby galaxies, in- 
cluding low-ionization, low-luminosity AGNs such as LINERs 
and transition objects, and even H II nuclei. The gas veloc- 
ity dispersions are comparable to, but typically somewhat less 
than, the stellar velocity dispersions measured over a phys- 
ical scale of typically 200-400 pc. The velocity discrep- 
ancy observed in H II nuclei and AGNs is about the same, 
((Jg/cr*) « 0.8, with a surprisingly small standard deviation of 
only - 0.2 dex. As in the work of Whittle (1992b, 1992c) and 
Nelson & Whittle (1996), we find that powerful radio jets, in the 
minority of nearby sources where they are found, do sometimes 
seem to exert secondary, super-virial acceleration on the NLR 
gas, although the correspondence is hardly one-to-one. How- 
ever, the main new insight obtained from our work is that the 
impact of the central AGN may be more widespread and sys- 
temic than previously suspected. We find that, not only among 
Seyfert galaxies but also in the characteristically less power- 
ful LINERs and transition nuclei, the offset between <r g and 
cr* varies mildly but systematically with the level of AGN ac- 
tivity. The effect can be seen with nuclear activity as gauged 
through optical (Ha) luminosity (Fig. 3a), radio luminosity 
(Fig. 4), and especially Eddington ratio (Fig. 5). The weak- 
est AGNs in the Palomar sample (LINERs and transition nuclei 
with L Ha w 3 x 10 37 erg s _1 , P bcm <; 10 18 W Hz" 1 , and L bol /L EdA 
w 3 x 10~ 6 ) are characterized by a g /a*~ 0.6-0.8, while lumi- 
nous sources (Seyferts with Lu a ~3x 10 41 erg s" 1 , P(, cm w 10 23 
W Hz" 1 , and L bo i/L E dd ~ 4 x 10" 2 ) can attain cr g /a^ 1.2-1.4. 
Remarkably, H II nuclei, even those spanning a similar range 
in Ha luminosity and Hubble type as the AGNs, do not fol- 
low these trends (Fig. 3b). Evidently the physical nature of the 
activity — stellar or nonstellar — matters. 

It is worth noting that Whittle (1992b, 1992c) and Nelson & 
Whittle (1996) find (a g /(T*) ~ 9 for their sample of Seyferts 
(after excluding the luminous, jetlike radio sources), somewhat 
larger than the average value of (cr g /cr*) w 0.8 for the Palomar 
objects. This difference is easy to explain in view of the lu- 
minosity dependence just mentioned, since Whittle's sample 
is more luminous than the Palomar objects. In their analy- 
sis of emission-line objects from the Sloan Digital Sky Sur- 
vey (SDSS), Chen et al. (2008) find, as we do in this study, 
that agio*, varies slightly but systematically with nuclear spec- 
tral classification: transition objects have the lowest values and 
Seyferts the highest, with LINERs in between. This trend re- 
flects none other than the dependence of a g /a* on luminosity 
or Eddington ratio. Ho (2008, 2009) proposes that the mass 
accretion rate, which scales with the Eddington ratio, is the pri- 
mary physical driver responsible for variations in spectral clas- 
sification. The accretion rate increases systematically along the 
sequence transition objects — > LINERs — > Seyferts. 

If the size of the NLR scales with luminosity, Laor (2003) 
suggests that the NLR in low-luminosity AGNs may be suffi- 
ciently compact that its kinematics may be influenced more by 
the gravitational potential of the central black hole than by that 
of the surrounding bulge. Spatially resolved observations of the 
NLR in some LINERs at Hubble Space Telescope resolution in- 
deed do find that the line widths are largest interior to the black 
hole's gravitational sphere of influence, gradually merging at 
larger radii with the stellar dispersion of the bulge (Walsh et al. 



2008). If this effect were dominant, however, u g la^ would rise 
with decreasing luminosity, exactly the opposite of what is seen 
in our sample. Evidently on a scale of several hundred parsecs 
the integrated kinematics of the NLR are controlled primarily 
by the gravitational potential of the stars and by AGN feedback. 

The discovery of a tight relation between central black hole 
mass and bulge stellar velocity dispersion (the M B h-c* relation: 
Gebhardt et al. 2000; Ferrarese & Merritt 2000) has prompted 
renewed interest in efficient methods to estimate 0%, especially 
in objects such as luminous AGNs where this parameter is chal- 
lenging to measure directly (Greene & Ho 2006). In view of 
the approximate equality between a g and 0% and the ease with 
which [O III] A5007 can be detected in AGNs, Nelson (2000) 
suggested that the stellar velocity dispersion of the host's bulge 
can be estimated by assuming cr* = cr g ([0 III]). Numerous stud- 
ies have since adopted this strategy to investigate the M Bli -&* 
relation in AGNs (e.g., Wang & Lu 2001; Grupe & Mathur 
2004) and its possible evolution with redshift (e.g., Shields et 
al. 2003; Salviander et al. 2007). However, the large observed 
scatter of the cr g -<7* correlation (Boroson 2003), not to mention 
the kinematic anomalies long known to afflict the [O III] line 
(Boroson 2005), gave cause for concern. To mitigate these un- 
certainties, Greene & Ho (2005a) recommended that only the 
core of the [O III] line be used, and, if possible, that [O III] 
be abandoned altogether in favor of lower-ionization transitions 
such as [O II] A3727 or [S II] AA6716, 6731. 

This study shows that the width of [N II], another low- 
ionization line, also traces the stellar velocity dispersion of the 
bulge. The scatter is modest, ~ 0.2 dex, but not insignificant. 
Since M BH oc cr* 4 (Tremaine et al. 2002), a 0.2 dex uncertainty 
on <t* translates into a 0.8 dex (or factor of ~6) uncertainty 
on Mbh- Importantly, we find that the velocity dispersion of 
the ionized gas does not identically match the velocity disper- 
sion of the stars, but rather, on average, o^/er^ss 0.8. Closer 
examination reveals that the residuals of the a g -a* relation in 
fact depend on AGN properties. The principal "second pa- 
rameter" appears to be either the Eddington ratio (Lboi/^Edd) or 
some measure of the AGN luminosity (we use Ln a ). The resid- 
uals seem to correlate somewhat better with Lboi/^Edd than with 
L Ha , although from our data alone we cannot truly tell which 
variable is more fundamental. Since the optical NLR luminos- 
ity correlates strongly with radio continuum luminosity (e.g., 
Ho & Peng 2001; Ulvestad & Ho 2001; Nagar et al. 2005), 
it is hardly surprising that Act correlates with radio power too, 
although the scatter is larger and the statistical significance is 
inferior compared to either Ln a or Lboi/^Edd- The dependence 
of the residuals of the a g -a* relation on Lboi/^Edd was first no- 
ticed by Greene & Ho (2005a; later confirmed by Bian et al. 
2006) for a large sample of more luminous AGNs selected from 
SDSS, and, within the uncertainties, the Palomar sources seem 
to follow roughly the low-luminosity extrapolation of the SDSS 
sample (Fig. 5). Curiously, Greene & Ho did not see any cor- 
relation between Ac and AGN optical ([O III]) luminosity or 
radio power. It is possible that the large linear scales probed 
by SDSS (the 3" fibers correspond to 5.4 kpc at z = 0.1) intro- 
duce substantial uncertainties into measurements of the bulge 
and NLR. 

The dependence of Act on AGN properties can be exploited 
to sharpen the er g -<7* correlation as a tool to predict er» when 
the latter is, as is often the case, difficult or impractical to mea- 
sure directly in AGNs. Even when a low-ionization line or only 
the core of the [O III] line is being used, as recommended by 
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Greene & Ho (2005a), it is important to realize that o g lo* is not 
identically equal to unity. Rather, o^lo* varies from ~ 0.6 to 
~ 1 .4 as a function of luminosity or Eddington ratio. Equations 
(1) and/or (3) should be used to obtain a more refined estimate 
of cr*. Moreover, significant scatter still remains (<~ 0.15 dex), 
even after correcting for the second parameter, most of which 
may be intrinsic and irreducible owing to the complexities of 
the NLR (Rice et al. 2006). 

4.2. Evidence for AGN Feedback 

From a physical point of view, what is the origin of the Aa- 
L Ha or Acr-Lboi/^Edd correlations? As discussed further below, 
if the gas derives principally from mass loss from bulge stars, 
its kinematics should generally track the kinematics of the stars. 
But because the gas is collisional and experiences hydrodynam- 
ical drag against the surrounding hot medium, we expect it to 
be kinematically slightly colder than the stars. In the absence of 
additional energy input from other sources, we anticipate a g /a* 
<; 1, as observed. As additional energy is injected into the sys- 
tem, for example from activation of the central black hole, the 
gas gains energy, to the point that a g approaches or even over- 
takes <7*. Precisely how this is accomplished is unclear. Do 
AGNs impart mostly mechanical energy from outflows/winds 
or more highly collimated fast jets? Does radiation pressure 
matter? The empirical trends presented in this paper offer some 
insights into these issues, which ultimately have bearing on cur- 
rent interests in AGN feedback in galaxies. 

The fact that the emission-line widths become systematically 
broader with increasing AGN activity — whatever the principal 
source of energy — implies that the gas is primarily accelerated 
rather than directly heated. The Ha luminosity, and presumably 
also the bolometric luminosity, strongly correlates with Act, but 
it seems doubtful that the radiation density itself plays a central 
role because H II nuclei with the same range in luminosity and 
sampled over a similar volume do not show the same effect. 
Compared to nuclear star formation, AGNs have a characteris- 
tically harder ionizing spectrum and a more centrally concen- 
trated radiation field. Whether these factors result in more effi- 
cient radiative acceleration of the ionized gas in AGNs needs to 
be further investigated. 

The radio-emitting plasma is an obvious suspect — that is, a 
scaled-down, but qualitatively similar agent as that responsible 
for super- virial acceleration in the most extreme, jet-dominated 
sources. This is a tenable hypothesis because nearly all Seyferts 
(e.g., Kukula et al. 1995; Ho & Ulvestad 2001; Ulvestad & Ho 
2001) and LINERs (Nagar et al. 2005) and a sizable fraction 
of transition objects (Filho et al. 2000, 2002) contain nuclear 
radio sources, many with morphological structures resembling 
jets. The radio luminosity function of low-luminosity AGNs 
varies continuously over at least 4 orders of magnitude in ra- 
dio power, from P 2Qcm « 10 19 to 10 23 W Hz" 1 (Ulvestad & 
Ho 2001; Filho et al. 2006). There is no characteristic power 
that might demarcate a luminosity threshold above which ra- 
dio sources become effective in perturbing the NLR. Although 
the degree of radio-loudness (expressed as the relative fraction 
of the radiative output emerging in the radio band) actually de- 
creases with increasing Eddington ratio (Ho 2002; Terashima 
& Wilson 2003; Greene et al. 2006), Ho (2008) notes that ex- 
tended, jetlike structures are more prevalent in high-Lboi/^Edd 
sources (Seyferts) than in low-Lboi/^Edd sources (LINERs). The 
radio cores in LINERs, despite being energetically more dom- 
inant, may be less effective at stirring up the NLR because 



of their compact structure. The extended, linear features in 
Seyferts, on the other hand, can impact a larger pool of the 
circumnuclear material. Jet-induced interactions have been in- 
voked to explain the morphology and kinematics of the NLR in 
a number of well-studied Seyfert galaxies (e.g., Mrk 78: Whit- 
tle & Wilson 2004; NGC 1068: Das et al. 2006; NGC 4151: 
Mundell et al. 2003). Despite the natural appeal of radio jets, 
however, our analysis shows that radio power correlates more 
poorly with Act than Ha luminosity or Eddington ratio. 

We speculate that the dominant source of energy input comes 
from accretion disk winds, a common feature in AGNs (Cren- 
shaw et al. 2003). Disk winds or outflows presumably agi- 
tate the NLR through shocks, provided that they do not violate 
other constraints that limit the role of shocks in the excitation 
of the gas (Ho 2008). The ability of an accretion disk to gen- 
erate a wind may depend critically on Lboi/^Edd (Proga 2007), 
and systems such as LINERs with ultra-low Lboi/^Edd seem to 
lack outflows altogether (Ho 2008). Consistent with Greene 
& Ho (2005a), we propose that the Eddington ratio primarily 
accounts for the secondary line broadening in the NLR. The 
apparent interchangeability between Lboi/^Edd and L Hq is prob- 
ably an artifact of the small dynamic range in black hole mass 
for the majority of the Palomar sample, coupled with the sizable 
intrinsic scatter in the Og-o* relation. 

4.3. Line Width-Luminosity Correlation in AGNs 

Phillips et al. (1983) first noticed that the luminosity of the 
[O III] line in Seyfert galaxies loosely correlates with its line 
width. This result has since been extended to much larger sam- 
ples of Seyferts by Whittle (1985, 1992b) and Gu et al. (2006), 
as well as to lower-luminosity AGNs in the Palomar survey by 
Ho et al. (2003; using Ha for line luminosity and [N II] for 
line width). The form of the correlation is L N lr °c FWHM", 
where, depending on the study and sample, the slope can take 
a wide range of values, from a m 3 to 6, with a typical value 
of <~ 4. The physical origin of the line width-luminosity rela- 
tion has been unclear. Whittle (1992b) suggests that the line 
width-luminosity relation is actually a secondary correlation, 
reflecting on the one hand the more primary link between bulge 
mass and NLR velocity, and on the other various possible inter- 
dependences among bulge mass, NLR luminosity, radio power, 
and line width. 

In light of the results of this study and other recent develop- 
ments, we offer a simple explanation. We have confirmed, as 
has long been suspected, that the widths of the narrow emis- 
sion lines largely reflect the gravitational potential of the bulge, 
such that FWHM oc cr* . The bulge, meanwhile, directly links to 
the mass of the central black hole via the fundamental relation 
Mbh oc a\ (Tremaine et al. 2002). Now, a given Mbh radiates 
up its Eddington luminosity, which in turn is related to the NLR 
luminosity by a bolometric correction and the Eddington ratio 
(accretion rate) of the system. Hence, AGNs should populate a 
distribution whose upper envelope follows a ridgeline defined 
by L NLR oc FWHM 4 . At a given FWHM, Lnlr has a maximum 
but no hard minimum (other than that dictated by detection sen- 
sitivity). Indeed, this is what is seen. Ho et al. (2003) presented 
the line width-luminosity correlation separately for each of the 
three classes of AGNs in the Palomar survey, and it is clear that 
the zero point of the correlation shifts with the average lumi- 
nosity of each class. 

4.4. Source of the Gas and Its Dynamical Support 
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As discussed in Ho (2008, 2009), the gas that ultimately fu- 
els AGNs in nearby galaxies can be readily supplied internally 
by mass loss from evolved stars (red giants and planetary neb- 
ulae) in the bulge. The same argument, the gist of which was 
already advanced by Minkowski & Osterbrock (1959) to ex- 
plain the origin of nebular emission seen in ellipticals, can be 
extended to include the entire gas reservoir that makes up the 
NLR in low-luminosity AGNs. Assuming T e = 10 4 K, n e = 100 
cm" 3 , Ha/H/3 = 3.1, and an effective recombination rate for Ha 
of af a = 9.36 x 10" 24 cm 3 s" 1 (Osterbrock 1989), 



M NLR = 2.97 x 10 3 



100 cm" 



Lhc, 



10 38 



ergs 



Mr. 



(4) 



From the population statistics of the Palomar survey assembled 
in Ho et al. (2003), LINERs and transition objects have me- 
dian Lh q ~ 7 x 10 38 erg s" 1 and n e ps 200 cm" 3 , and hence 
A/nlr ~ 1 x 10 4 Mq. The corresponding values for Seyferts are 
L Ha w 40 x 10 38 erg s" 1 , n e « 400 cm" 3 , and M NL r « 3 x 10 4 
Mq. Strictly speaking, these values only pertain to the central 
200 pc x 400 pc region sampled by the spectroscopic aper- 
ture, but given the high degree of central concentration of the 
line-emitting gas (Ho 2008), they probably underestimate the 
integrated values for the entire NLR by no more than a fac- 
tor of ~ 2. These mass estimates can be compared with the 
amount of material shed from evolved stars. For a Salpeter stel- 
lar initial mass function with a lower-mass cutoff of 0. 1 Mq, an 
upper-mass cutoff of 100 Mq, solar metallicities, and an age of 
15 Gyr (Padovani & Matteucci 1993), 



M, w3x 10" 



M yr 



(5) 



With a median Hubble type of Sa, the host galaxies of the Palo- 
mar AGNs have sizable bulges and correspondingly healthy gas 
supply rates. For a median bulge luminosity of M B w -19.3 
mag (Ho et al. 2003) and B-V = 0.8 mag (Fukugita et al. 
1995), M.J, m 0.5M Q yr -1 . These values pertain to the entire 
bulge. Within the spectroscopic aperture the nuclear stellar con- 
tinuum luminosities (Ho et al. 1997a) are approximately a fac- 
tor of 10 lower than the integrated luminosities, which implies 
rj 0.05 M Q yr . Thus, stellar mass loss can sustain the gas 
reservoir in the NLR if the stellar debris survives for periods 
longer than ~ 2 x 10 5 - 10 6 yr before it dissipates and merges 
with the surrounding hot interstellar medium (Mathews 1990; 
Parriott & Bregman 2008). 

If the ionized gas is mainly produced internally through stel- 
lar mass loss in the bulgcl what is the expected kinemati- 
cal signature of this material? The supersonic motion of the 
ejecta through the ambient medium generates a strong shock 
that thermalizes the gas to the kinetic temperature of the stars, 
- 10 6 - 10 7 K, on a timescale of - 10 s - 10 6 yr (Sanders 1981; 
Mathews 1990). Not all of the gas gets heated, however. Nu- 
merical simulations of mass loss from evolved stars in ellipti- 
cal galaxies, after accounting for radiative cooling, predict that 
some fraction of the gas, perhaps ~ 20%, remains cool (Parriott 
& Bregman 2008). Photoionization of this residual gas, as well 

4 We do not mean to suggest that external acquisition of gas does not take place. Indeed, an external origin has been invoked to explain the distribution of kinematic 
axis misalignments between the gas and stars (Bertola et al. 1992; Caon et al. 2000; Sarzi et al. 2006), as well as certain H I properties (Ho 2007), in early-type 
galaxies. 

5 The optical line emission in ellipticals and bulges in principle can be excited by a variety of mechanisms, including photoionization by young or old hot stars, 
photoionization by a central AGN, shocks, self-irradiation from cooling condensations, and even cosmic ray heating by the radio-emitting plasma. With the possible 
exception of the weakest emission-line objects, however, Ho (2008) argues that the majority of the low-luminosity AGNs in the Palomar survey — essentially all 
nearby galaxies with sizable bulges — are powered by nonstellar photoionization. 



as the fresh ejecta prior to thermalization, produces the nebular 
emission visible in the optical|f| During the early stages of their 
evolution, the mass loss envelopes orbit roughly cospatially 
with their parent stars and thus closely track their natal veloci- 
ties. As bulge stars do not have a strong rotational component 
to their velocities, neither should the gas shed by them. After 
normalizing the gas line widths by the observed (projected) ro- 
tational velocity of the large-scale (H I) disk (Fig. 6), the signif- 
icant, positive trend with galaxy inclination angle that remains 
implies that the ionized gas does not reside in a plane aligned 
with the large-scale galaxy. Moreover, the ejecta quickly lose 
memory of their birth sites as the material experiences drag de- 
celeration against the surrounding hot medium. As discussed 
by Mathews (1990), the warm, optical-emitting clouds cannot 
be self-gravitating. Instead, they must be pressure-confined 
by the ambient hot medium, and hence their kinematics will 
largely imprint the local velocity of the hot gas. Since the hot 
gas is in virial equilibrium with the bulge stars, to first order 
we anticipate a g sa ov In detail, our analysis finds a g w 0.8cr*. 
As the referee points out, supersonic turbulence in the hot gas 
rapidly dissipates through shocks. It is thus natural for a g to 
be slightly subvirial because the turbulence — the source of the 
random motions — is subsonic. When an AGN is present, the 
hot gas may be more agitated. 

The scenario presented above provides a useful framework 
for interpreting a number of trends gleaned from spatially re- 
solved observations of the central regions of galaxies. Hubble 
Space Telescope images frequently reveal dust lanes, spiral pat- 
terns, and a variety of other small-scale fine features in galaxy 
centers, but rarely do the structures ever take the form of thin, 
well-defined disks (Phillips et al. 1996; Carollo et al. 1997; 
Malkan et al. 1998; Pogge et al. 2000; Ho et al. 2002; Martini 
et al. 2003; Simoes Lopes et al. 2007). Since the dust generally 
traces the nebular gas (e.g., Ravindranath et al. 2001; Falcon- 
Barroso et al. 2006), we can infer that the ionized gas more or 
less follows the same topology as the dust. The complex mor- 
phology of the gas strongly hints that its velocity field must be 
similarly chaotic. Where available, observations support this 
picture. For example, a number of authors have noted that the 
central rotation curves of the ionized gas often rise more slowly 
than the circular velocities predicted from the luminosity pro- 
file of the stars (Caldwell 1984; Caldwell et al. 1986; Fillmore 
et al. 1986; Kent 1988; Kormendy & Westpfahl 1989; Cinzano 
& van der Marel 1994; Bertola et al. 1995; Fisher 1997; Cin- 
zano et al. 1999; Pignatelli et al. 2001), suggesting that random 
motions contribute appreciably to the dynamical support of the 
gas. On scales probed by the Hubble Space Telescope, very few 
objects exhibit clean signatures of dynamically cold disks un- 
dergoing circular rotation (Sarzi et al. 2001; Ho et al. 2002; 
Atkinson et al. 2005; Noel-Storr et al. 2007). Although the 
statistics are still quite limited, evidence is emerging that the 
degree of kinematic disturbance may be directly linked to the 
level of AGN activity (Verdoes Kleijn et al. 2006; Dumas et al. 
2007; Walsh et al. 2008; § 4.2). 

We end with the following note. The above discussion pre- 
supposes that the optical line emission comes exclusively from 
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stellar ejecta prior to their being mixed into the hot phase. 
There may be another component. Once the ejecta become 
thermalized, the hot gas, at least in the inner cores of bulges, 
should cool and condense on a relatively short timescale. Re- 
cent Chandra and XMM-Newton observations of ellipticals and 
bulges show that the diffuse X-ray-emitting medium typically 
has temperatures of kT « 0.3- 1 keV and densities of n « 
0.1-0.3 cm" 3 (see Ho 2009 and references therein). Assum- 
ing kT = 0.5 keV and n = 0.1 cm" 3 , the cooling time is merely 
~6x 10 7 yr. This simple picture of galactic-scale cooling flows 
encounters a number of long-standing difficulties familiar in the 
context of galaxy clusters, and their resolution appeals to sim- 
ilar remedies (Mathews & Brighenti 2003). A key ingredient 
to counter the overcooling problem invokes some form of en- 
ergy feedback. It is tempting to associate the evidence for AGN 
feedback discussed in this paper (§ 4.2) with such a mecha- 
nism. The kinematics that ensue from the cooling flow scenario 
are likely to be quite complicated, but in general we also ex- 
pect the cooling condensations to be in pressure balance with 
the surrounding hot medium and hence to reflect its kinematics. 
Depending on the shape of the potential, the cooling flow fila- 
ments can settle in different planes at different radii (Tohline et 
al. 1982). This may account for the variety of dust structures 
often seen in the central regions of early-type galaxies (e.g., 
Lauer et al. 2005). 

5. SUMMARY 

We use a new catalog of stellar velocity dispersions to study 
the relationship between the kinematics of the ionized gas and 
the stars in the central few hundred parsecs of a large sample of 
nearby galaxies. The gas dispersions are based on the width of 
the [N II] A6583 line. Robust measurements of both a g and cr* 
are available for 345 galaxies spanning a wide range in Hubble 
type and level of nuclear activity. The principal results of our 
analysis can be summarized as follows. 

1 . The velocity dispersion of the ionized gas strongly cor- 
relates with, and is roughly comparable to, the velocity 
dispersion of the stars over the velocity range ~ 30-350 
km s" 1 , across all Hubble types. This confirms the no- 
tion that the widths of the narrow emission lines primar- 
ily reflect the virial velocity of the gravitational potential 
of the stellar bulge. In detail, a g /a* ranges from <~ 0.6 
to 1.4, with an average value of 0.80. 

2. The ratio cr g /(T* shows no strong dependence on Hubble 
type, for galaxies ranging from E to Sbc. Among galax- 
ies of type Sc and later, systems with < 40 km s" 1 , 



all spectroscopically classified as H II nuclei, a g la^ ;> 
1 . The ionized gas maintains a minimum "floor" disper- 
sion of o g i=a 30 km s" 1 , which might signify a lower 
threshold level of interstellar turbulence generated by 
energy feedback from massive stars. The o g -o* relation 
does not depend on the presence of a large-scale bar or 
on the local environment of the galaxy. 

3. The residuals of the o g -o* relation for AGNs (Seyferts, 
LINERs, and transition objects) correlate strongly with 
the nuclear optical (Ha) luminosity, and perhaps even 
stronger with the Eddington ratio. Radio power plays 
a less essential role. We argue that energy injected by 
accretion disk winds or outflows, not radio jets, provide 
a source of AGN feedback that accelerates the NLR gas 
to speeds above the virial velocities of the stars. 

4. We confirm that the line widths of narrow emission lines 
can be used to obtain reasonably reliable estimates of 
cr*, but, for maximum accuracy, preference should be 
given to low-ionization transitions such as [N II] and 
care should be taken to correct for the secondary de- 
pendence of a g on AGN luminosity or Eddington ratio. 

5. We show that the mass budget of the ionized gas can 
be accounted for by mass loss from evolved stars. The 
optical line emission arises from photoionization of the 
just-shed, mass loss envelopes prior to their being ther- 
malized, in conjunction with an additional contribution 
from filaments recondensing from the hot medium. We 
argue that the line widths are dominated by random mo- 
tions, reflecting the kinematics of the hot gas that pres- 
sure confines the line-emitting clouds. 

6. Lastly, we suggest that the line width-luminosity corre- 
lation in AGNs actually defines the upper envelope of 
a distribution of points governed by the <jg-o* relation, 
the Mbh-c* relation, and the Eddington limit. 
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